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I briefly review advances in the understanding and modeling of relativistic stellar dynamics 
around massive black holes (MBHs) in galactic nuclei, following the inclusion of coherent 
relaxation and of secular processes in a new formal analytic description of the dynamics. 


1 Relaxation in galactic nuclei 

The dense, centrally concentrated stellar cluster that exists around most MBHs offers oppor¬ 
tunities for strong, possibly destructive interactions between it and stars. These include direct 
plunges, leading to tidal disruption flares or gravitational waves (GW) flares, inspiral processes 
leading to quasi-periodic GW emission from extreme mass ratio inspiral events (EMRIs) 1 , tidally 
powered stars (“squeezars”) 2 , strong tidal scattering 3 or capture by massive accretion disks. 
These processes affect MBH growth and may create exotic stellar populations around MBHs 3 . 

This naturally leads to the question “How do stars closely interact with, and fall into a 
MBH, and at what rates?” This is known as the stellar dynamical “loss-cone problem”. It is a 
non-trivial problem, in spite of the presence of so many stars so close to the MBH, because the 
phase space volume of unstable orbits is minute. The few stars initially on such orbits quickly 
fall into the MBH on the short dynamical timescale, and then the rates would drop to zero, 
if it were not for dynamical processes that deflect additional stars from stable orbits to those 
with velocity vectors that point toward the MBH, within the loss-cone (Fig. 1 left). Thus, the 
loss-cone question is essentially the question: “how do galactic nuclei randomize and relax?” 

1.1 Non-coherent 2-body relaxation (NR) 

The discreteness of stellar systems leads to non-coherent 2-body relaxation (NR) (Fig. 1 left). 
This guarantees a minimal relaxation rate, on a timescale Tnr Q 2 P(r)/N*(r) log Q, where 
Q = M,/M* is the MBH/star mass ratio, P(r) is the radial orbital period, and N+(r) the number 
of stars inside r. Because the impact parameter b of these point-point interactions can be small, 
NR is boosted by the Coulomb factor log(6 max /6 m in) = log Q. Trtr is the timescale for changes 
of order unity in energy, Tr. It is however easier to drive a star into the MBH by reducing its 
angular momentum L and making its orbit more radial, than by reducing the orbital energy 
E < 0, and shrinking the orbit. The timescale for changing j = LjL c (a ) = \f \ — e 2 from j to 0 
is T l = j 2 T E (L c = JGM^l is the circular L, a the sma and e the eccentricity). 

In the absence of dissipation, stars with j < 1 are deflected by L-scattering at nearly constant 
a to the innermost stable orbit (ISO), at ji so = 4 \JryJa ( r g = GM./c 2 ) and then plunge directly 
into the MBH 5 (Fig. 1 center). When a dissipative mechanism is present (e.g. GW), phase 
space is divided in two (Fig. 1 left). Below some critical sma a c , all stars eventually cross the 




Figure 1 - The 2-body relaxation-driven loss-cone. Left: A star is scattered to an orbit in the loss-cone, which 
takes it close enough to the MBH for a strong (possibly destructive) interaction with it. Center: The loss-cone 
phase space in terms of the normalized angular momentum j = L/L c = \/l — e 2 and the semi-major axis a, 
without dissipation. Right: The same, but with a dissipative mechanism (here, the emission of GWs). See text. 


“inspiral line” where id-dissipation is faster than j-scattering, and then inspiral gradually into 
the MBH as EMRIs. Stars above a c plunge directly. The respective rates of plunges and inspirals 
can then by estimated by the ratio of number of stars on the relevant scales (the MBH radius 
of influence r / 1 for plunges (e.g tidal disruptions), and a c for inspirals) over T^ R on that scale: 
R p ~ N+(r h )/T NR \og(L c /L iso ) and ~ N+(a c )/T NR log(L c /L iso ). Because iV*(a c ) < N^rh), 
the inspiral rate is much lower than the plunge rate 6 , typically Ri ~ 0.01i? p . 

1.2 Coherent resonant relaxation (RR) in nearly-spherical systems 

Resonant relaxation 4 is a process of rapid L-relaxation that occurs when the gravitational 
potential is symmetric enough to restrict the evolution of orbits on timescales much longer than 
the orbital time (e.g. nearly-fixed Keplerian ellipses in the nearly-Keplerian potential close to 
a MBH, where the stellar mass is negligible, but far enough so that GR effects are weak). In 
that case, a test orbit will feel a residual torque from the static, orbit-averaged background of 
stellar “mass wires”, which persists for a long coherence time T c , until small deviations from 
symmetry accumulate and randomize the background. These orbit-orbit interactions randomize 
the angular momentum on a timescale T RR rs j Q 2 P(r) 2 /N±(r)T c . Unlike NR, these extended 
objects do not undergo close interactions. Rather, RR is boosted by the long coherence time. 

RR is relevant for the loss-cone problem because it is possible to have T c P in the 
symmetric potential near a MBH, so that T rr /Tn r ~ (log Q)P/T C <C 1. That is, angular 
momentum evolution, and in particular that leading to j —)■ 0 and strong interactions with the 
MBH, can be greatly accelerated. Unchecked, RR will completely suppress EMRIs by driving 
all stars into plunge orbits (Fig. 2 center). However, very eccentric orbits undergo GR in¬ 
plane (Schwarzschild) precession, which quenches RR by rapidly alternating the direction of 
the residual torque on the orbit. This motivated the “fortunate coincidence conjecture” 7 : The 
0(/3 2 j~ 2 ) GR precession becomes significant before 0(f3 5 j~ 7 Q~ 1 ) GW dissipation, and this may 
allow EMRIs to proceed unperturbed, decoupled from the background stars. 

1.3 The Schwarzschild Barrier 

The first full PN2.5 )V-body simulations 8 revealed a surprising result: not only does GR pre¬ 
cession quench RR before the GW-dominated regime, as conjectured, but there appears to be 
some kind of barrier in phase space, dubbed the Schwarzschild Barrier (SB), which prevents 
the orbits from evolving to j —> 0. Instead, they appear to linger for roughly T c near the SB, 
where their orbital parameters oscillate at the GR precession frequency, and then they evolve 
back to j —> 1. An early analysis 9 suggested that this behavior is related to precession under 
the influence of a residual dipole-like residual force. However, a full self-consistent explanation 
of the SB was lacking, and its very existence and nature remained controversial. 













Figure 2 - The loss-cone phase space. Left: A schematic for a model of the Milky Way nucleus (M. = 4 x 10 6 M q, 
M* = 10 Mq). Orbits are unstable in the gray region left of the ISO line. Dynamics are dominated by GW 
dissipation inside the GW curve at the bottom right. The horizontal line tangent to the upper tip of the GW 
line is the critical sma separating plunge and inspiral tracks (one example shown for each). RR dominates over 
NR in the shaded region in the middle right. RR is ineffective on short timescales below the diagonal AI line 
just below the RR region. Center: The phase space density and loss rates calculated by MC simulations, for an 
artificial model without GR precession to quench RR, resulting in rapid plunges (circles at ISO line) and complete 
suppression of EMRIs. Right: The same, but with GR precession, which suppresses RR (RR remains strong 
only inside the Trr = 0.1, ITjvh contours). This enables EMRIs (triangles at bottom of GW region). See text. 


I now describe briefly a new formal framework for expressing coherent relaxation and secular 
processes in galactic nuclei 10 , and discuss implications for the steady state phase space structure 
and loss rates 11 . 


2 Hamiltonian dynamics with correlated background noise 


Two key insights inform the new advances in understanding coherent relativistic dynamics, which 
lie in the difhcult-to-treat interface between deterministic Hamiltonian dynamics and stochastic 
kinetic theory. (1) The effect of the background stars on a test stars should be described by a 
correlated noise model rf(t), whose degree of smoothness (differentiability) determines dynamics 
on short timescales 10 . (2) The long-term steady-state remains (unavoidably) the maximal 

entropy configuration, irrespective of the details of the nature of the relaxation processes 11 . 

This formal treatment of PN1 dynamics in the presence of correlated (RR) noise rj (a 3- 
vector in L-space) allows to write a phase-averaged leading-order (£ = 1) Hamiltonian j-L\ and 
derive stochastic EOMs for the orbital elements of a test star, x = (j, </>, cos 9) and the argument 
of periapse if, which precesses at frequency v p (j), 

X = v T , x (x,if) • r? , ip = v T ^(x,if)-ri + v p (j), (1) 


where v T is the RR torque frequency. This r] formalism allows to evolve a test star in time for a 
given realization of the noise. Moreover, even though r/ is time-correlated, it is possible to derive 
(and validate with the stochastic EOMs) approximate diffusion coefficients (DCs) 2 , which 
allow to evolve in time the probability density p(j) with the Fokker-Planck (FP) equation, 
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J~C(t) is the Fourier transform of rj's auto-correlation function (ACF). The explicit dependence 
of D 2 on the spectral power of the noise at the precession frequency is an expression of adiabatic 
invariance (AI). If, and only if the noise has an upper frequency cutoff, as it must if it is smooth 
(this is physically expected, since the background noise is the superposition of continuous orbital 
motions), then there is a critical jo such that for j < jo the precession is fast enough so that 
D 2 (j) —> 0, and the star decouples from the background resonant torques (Fig. 3 left, center). 
This describes the dynamics of the SB: it is not a reflecting boundary, but a locus in phase space 
where diffusion rapidly drops due to AI. Since diffusion to yet lower j slows further down, while 
diffusion to higher j speeds further up, orbits statistically seem to bounce away from the SB. 












































































































Figure 3 - The smoothness of the noise model and diffusion dynamics. Left: Three r] models and their ACF: 
discontinuous steps (C 0 ,), continuous but not continuously differentiable (C 1 with exponential ACF), smooth 
(C°° with Gaussian ACF). Center: The corresponding Do; note the steep cutoff at j < 0.1 for the smooth noise 
model. Right: The MC simulations of j-only evolution reproduce the AI/SB limit at j < jo in the absence of 
NR, but NR erases this feature completely on timescale t —» Tnr- 


3 The steady state loss-cone 

NR is impervious to AI. When t —> T/vr, the system approaches the maximal entropy solution 
( dN/dj = 2 j, Fig. 3 right). Monte Carlo (MC) simulations of the probability density, branching 
ratios and loss rates with the ij formalism (Fig. 2 right) show that the RR-dominated region in 
phase space is well separated from the plunge and inspiral loss-lines, so the effect of RR on the 
loss rates is small (< x2 — 3). Specifically, we conclude that GR quenching of RR is effective, 
so the EMRI rates remain largely unaffected by RR. RR can be significant for processes whose 
loss-line crosses the RR-dominated region, e.g. destruction by interaction with an accretion disk. 

4 summary 

NR, RR, GW dissipation and secular precession can be treated analytically as effective diffusion 
with correlated noise. The steady state depends mostly on NR, which erases AI. RR can be 
important in special cases. The rj formalism provides stochastic EOMs for evolving test particles 
and an FP/MC diffusion procedure for evolving the probability density. This makes it possible 
to model the relativistic loss-cone in galactic nuclei with realistic IV* 1, (unlike direct IV-body 
simulations), and obtain the branching ratios, loss rates and steady state stellar distributions. 
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